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The structure of the Cooper minima in the transition probabilities and photoionization cross- 
sections for low-excited and Rydberg nS, nP, nD and nF states of alkali-metal atoms has been 
studied using a Coulomb approximation and a quasiclassical model. The range of applicability of the 
quasiclassical model has been defined from comparison with available experimental and theoretical 
data on dipole moments, oscillator strengths, and photoionization cross-sections. A new Cooper 
minimum for transitions between rubidium Rydberg states has been found. 
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I. INTRODUCTION 

Spectral line series of alkali metal atoms display re- 
markable features with prominent minima in the tran- 
sition probabilities, emission oscillator strengths or pho- 
toionization cross-sections [H4l3|. These minima arise 
from the cancellation of the radial integral for some 
transitions, depending on the overlap between the wave- 
functions of the initial and final quantum states of the 
atoms [j], [|| , and are well known as Cooper minima [!, H[ ■ 

Observation of the Cooper minima in the photoion- 
ization cross-sections and in the transition probabilities 
of the discrete spectrum provides valuable information 
about the electronic structure of the atoms. Minima 
in photoionization cross sections were first found exper- 
imentally in Ref. [|[ and explained 20 years later by 
Cooper [3|. The sharp minima in the emission probabil- 
ities for some Rydberg states of alkali-metal atoms were 
first discussed by Theodosiou @. 

The experimental investigation of Cooper minima can 
be used for the verification of theoretical calculations of 
spectroscopic properties of atoms and molecules. For ex- 
ample, experimentally measured photoionization cross- 
sections for sodium ground state showing the Cooper 
minimum Q, were in good agreement with the theo- 
retical calculations of Aymar [8|, which confirmed the 
accuracy of the theoretical model. 

Cooper minima in the discrete spectrum are revealed 
as a suppression of the two-photon photoionization [10j | 
or sharp decrease of the emission oscillator strengths Q . 
The map of these minima could be valuable for system- 
atic studies of the processes which involve a large number 
of transitions, such as calculation of lifetimes of Rydberg 
atoms [lij , blackbody-radiation- induced photoionization 
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rates [TH, fl(| or collisional ionization cross-sections of 
cold atoms 



Photoionization of alkali-metal atoms recently at- 
tracted a lot of interest, as it has taken a central 
role in experiments with cold atoms in far-off-resonance 
traps [l8|, EH, in photoionization spectroscopy [2fJrl22l ] . 
measurement of oscillator strengths |2l| - [23| , photoioniza- 
tion cross-sections and lifetimes of excited atoms [3 [2ll — 
[29j . and photoionization of the Bose-Einstein conden- 
sate Hla3. 



Rydberg atoms with large principal quantum num- 
bers ?i~50-100 recently received attention due to the 
progress achieved in experiments with cold atomic sam- 
ples. These samples are often prepared in optical traps 
whose intense laser field could make the ionization life- 
times of Rydberg atoms extremely short. However, for 
certain wavelengths, ionization rates could be signifi- 
cantly reduced due to Cooper minima in the photoioniza- 
tion cross-sections Q . Therefore it would be most useful 
to exploit trapping schemes operated at wavelen gths dis- 
playing Cooper minima to avoid photoionization |T8l. [32j . 



In this paper we have calculated the radial matrix el- 
ements of arbitrary bound-bound, bound-free and free- 
free transitions between S, P, D and F states of alkali- 
metal atoms using the quasiclassical model by Dyachkov 
and Pankratov (DP model) [HI, [34|. In section II we 
present examples of the Cooper minima for bound-free 
and bound-bound transitions and discuss the accuracy of 
the theoretical method. Numerical results are presented 
in Section III as density plots, revealing the Cooper min- 
ima in bound-bound, bound-free and free-free transi- 
tions. 
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II. THE QUASICLASSICAL MODEL AND ITS 
APPLICABILITY 



Radial matrix elements of the electric dipole transi- 
tions between arbitrary atomic states (e.g., bound-bound 
or bound-free transitions) are required to calculate the 
spectroscopic properties of atoms, including oscillator 
strengths, lifetimes, photoionization cross-sections, and 
rates of collisional ionization. 

Although alkali-metal atoms have a single valence elec- 
tron, only states with small quantum defects exhibit truly 
hydrogen-like behavior. Due to non-hydrogenic charac- 
ter of alkali-metal atoms, the calculation of radial ma- 
trix elements remains a challenging task, since no exact 
analytical solution for arbitrary transitions is available 
yet [35|. The oscillator strengths for alkali- metal atoms 
can strongly deviate from the values for hydrogen. Accu- 
rate calculation of the radial integral for transitions be- 
tween states with small angular momentum are difficult 
because of the need to take into account the interaction 
of the valence electron with the atomic core. 

A method based on the Coulomb approximation re- 
lies on the idea that the Rydberg electron is localized 
mostly outside the atomic core, where the potential is 
Coulombic. In the Numeric Coulomb Approximation [361 ] 
the radial wavefunctions are obtained by solution of the 
Shrodinger equation with the exact energies of the alkali- 
metal quantum states, expressed through the quantum 
defect (Rydberg-Ritz formula, atomic units are used in 
this paper): 



1 



2n 2 eff 



(1) 



Here n e / / = n — /il is an effective quantum number, 
is the quantum defect. Quantum defect accounts for the 
penetration of the valence electron into the ionic core 
of a Rydberg atom. The quantum defects are used as 
input parameter for the calculations and the integration 
is truncated at the inner core radius. 

Alternative forms of the Coulomb approximation were 
developed in [37} • The Modified Coulomb Approximation 
(MCA) is a generalization of the analytical expression 
for the hydrogen radial integral for non-integer quantum 
numbers. It allows direct calculation of the radial matrix 
elements without numeric integration. 

Further simplification of the calculations in the 
Coulomb approximation is achieved by extension of the 
quasiclassical approximation to the states with low prin- 
cipal quantum numbers [33l . l3~i| . The radial matrix el- 
ements are expressed through transcendental functions, 
which help to avoid inaccurateness of the direct numerical 
integration. This method can substantially improve both 
the speed and reliability of the calculations. However, the 
validity of most quasiclassical models was restricted by 
transitions between neighboring excited states [38l - l42l |. 

In our previous works |16| we used the quasiclassi- 
cal model developed by Dyachkov and Pankratov [HI, [34| . 
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FIG. 1. (Color online). The calculated radial matrix el- 
ements for bound-bound (a) 60F — > nD transitions in ru- 
bidium; (b) 60F — ¥ nD transitions in cesium; (c) 60D — > nP 
transitions in potassium; (d) 60P — » nS transitions in lithium. 



Their original approach provides more precise values for 
the wave functions of the Rydberg and continuum states, 
compared to the other quasiclassical models. Good 
agreement with numeric results based on NCA model [36[ 
or various model potential methods [43|, [44[ is observed. 

Radial matrix elements for transitions between excited 
states of alkali-metal atoms display numerous features in 
their depenencies on n. For example, minima were re- 
vealed in the transition probabilities for nD — > n'F in 
cesium [6] and for nD — > n'P in potassium. Hoogen- 
raad et al. [Toj observed theoretically and experimen- 
tally a Cooper minimum for nS — > n'P transitions in 
lithium. To study the validity of the quasiclassical model 
for bound-bound transitions, we have calculated the tran- 
sition probabilities for nL n' L' transitions in alkali- 
metal atoms with n < 100 and L, L 1 — 1,2,3 using 
the DP quasiclassical model, which is depicted in Ap- 
pendix B. Figure [1] shows the dependencies of the ra- 
dial matrix elements on the principal quantum number 
n for 60-F — > nD transitions in rubidium and cesium 
[Fig- da) and Fig.QTb), respectively], 60D — > nP transi- 
tions in potassium [Fig. [He)] an< ^ 60P —> nS transitions 
in lithium [Fig. QJd)] . The observed minima for lithium, 
potassium and cesium are in agreement with the results 
of Refs. @, [ll|. The minima for nF — > n'D transitions 
in rubidium appear only for high n, and have not been 
located yet, to the best of our knowledge. These minima 
lie in the microwave region of about 150 GHz and could 
be studied using microwave spectroscopy [45| . 

Direct measurement of the radial matrix elements is 
of great importance for verification of the theory. How- 
ever, due to lack of available experimental data for tran- 
sitions between excited states of alkali-metal atoms, new 
measurements are required. In order to benchmark the 
model, we have earlier measured the reduced matrix ele- 
ment for the diffuse series of rubidium [23T ]. 

We observed the Autler Townes splitting in a sam- 
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FIG. 3. (Color online). Lithium 2S photoionization cross- 
section. Comparison of the quasiclassical calculations with 
the experiment 0] and theory [I?], 5H . 



FIG. 2. (Color online), (a), (c) Comparison of the calculated 
reduced dipole moments for (a) rubidium atoms with experi- 
ment [23,[2^] and (c) cesium atoms with experiment [2ll. |24| . 
(b), (d) Cooper minima in photoionization cross-sections of 
nS (b) rubidium and (d) sodium atoms. Solid curves - this 
work. Broken curves - theoretical calculations from Ref. S|: 



pie of ultra-cold Rb atoms using a 3-level ladder system. 
Briefly, we monitored the absorption of a weak probe 



laser scanned over the 55i/2 — > 5P3/2 whilst simulta- 
neously a strong coupling laser, locked to the 5P 3 / 2 — > 
nD 5 /2 transition illuminated the atoms. The strong cou- 
pling laser gave rise to two absorption peaks separated 
by the Rabi frequency of the atom-coupling laser interac- 
tion. Knowledge of the laser intensity allowed us to mea- 
sure the dipole matrix elements of the 5P 3 / 2 - ► nD 5/2 
transitions to within 7% accuracy. 

We have also compared the reduced dipole moments 
calculated using DP model with other available exper- 
imental data on diffuse series of rubidium [22], |23[ and 
cesium E l , [2 4l|. Good agreement with experiments of 
Refs. |22L 23| is confirmed in Fig. [2fa) for rubidium 
5P — > nD transitions. For cesium the theoretical values 
in Fig. [2fc) are in agreement with the experiment only 
for the lowest nD states and differ from the experimental 
values for higher n by a factor of two. Our theoretical 
results for cesium are, however, in excellent agreement 
with the previous calculations of Ref. [24| • The observed 
disagreement between experiment and theory can result 
from improper account for core polarization for heavy 
cesium atoms [46[. However, we expect that for transi- 
tions between the states with larger principal quantum 
numbers the accuracy of the semiclassical approximation 
will be significantly increased due to smaller interaction 
of the Rydberg electron with the atomic core. 

Calculated Cooper minima in the photoionization 
cross-sections of the nS states of rubidium and sodium 
are shown in Fig. [^b) and (d), respectively. The solid 
curves in Fig. [2fb) and (d) represent the photoionization 
cross-sections calculated using our model based on the 
DP method. In the case of sodium, these are compared 
with the quantum-mechanical calculations of Ref. Q, 



shown as broken curves in Fig. ^d). Good agreement 
is observed for sodium nS states with n < 5. For higher 
states the positions of the Cooper minimum are signif- 
icantly shifted. In Ref. [49| it has been argued that 
the discrepancy in the quasiclassical calculations could 
be corrected by adjusting the phase in the radial integral 
in order to compensate the phase shift of the radial wave- 
function from the value given by the quantum defect. At 
the same time, the reliability of the quasiclassical approx- 
imation is also expected to improve with the increase of 
the principal quantum number, and only experimental 
data could confirm the validity of the theory. 

In figure [3] we have compared the calculated photoion- 
ization cross-section for lithium 2S state with the exper- 
iment of Ref. and theory of Refs. [Hill. It is seen 
that our approach provides a better agreement with the 
experiment. 

We conclude that the DP model is suitable for cal- 
culation of the radial matrix elements and photoioniza- 
tion cross-sections with an accuracy better than a fac- 
tor of two for low n and much enhanced for higher ex- 
cited states, as confirmed by the good agreement be- 
tween experiment and theory in recent lifetime measure- 
ments dlllii. 



III. MAPS OF THE COOPER MINIMA 



Using the quasiclassical model of Dclone et al. [42|]i m 
Ref. [Tl| it has been shown that the radial matrix ele- 
ments for bound-bound, bound-free and free-free tran- 
sitions could be expressed in a universal way through 
the numerically calculated relative matrix elements R re i , 
multiplied by the appropriate normalization factors: 
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(2) 



The prefactor (4/3) 1/3 /V (1/3) = 0.4108 results from 
the asymptotic expression for the quasiclassical matrix 
elements for n — > n + 1 transitions [111 ]. Relative ma- 
trix elements R re i (EL —> E'L') introduced in Eq.Q are 
convenient as these are slowly varying functions of E and 
E' . The dependence of R re i (EL — > £"i') on the energy 
E' of the final state passes smoothly through the ioniza- 
tion threshold 1JLJ . The asymptotic \E' — E\ 1 depen- 
dence of the radial matrix elements is incorrect for transi- 
tions between neighboring states with E ss E' , where the 
dipole matrix elements rapidly increase [Tlj ] . In this case 
the radial matrix elements can be calculated numerically 
using a DP model 0, H, or NCA (3<|. 

We have calculated the relative matrix ele- 
ments R re i (E n L — > E n 'L'), Rrel (E n L -> E'L'), 
Rrei (EL — > E'L) for transitions between S, P, D 
and F states of alkali-metal atoms, starting from the 
ground state. The energies of the continuum states were 
taken within < E < 0.5 (atomic units) to extend the 
results of Ref. |11| to the area where Cooper minima are 
expected for alkali-metal atoms. 



A. Rubidium 

Rubidium and cesium atoms are widely employed 
in laser cooling experiments, and we shall discuss 
them in more detail. The relative matrix elements 
Rrei (Si/2 ~^ Pi/ 2) f° r rubidium are shown in Fig. [4] Fol- 
lowing Ref. [ill ] , we present our numerical results as den- 
sity plots. We use both ^-scaled and n-scaled plots, since 
the latter are more appropriate to the relative matrix el- 
ements for bound-bound and bound-free transitions from 
states with large principal quantum numbers n. The 
signs of the radial matrix elements are unimportant in 
the calculation of transition probabilities and photoion- 
ization cross-sections, therefore we present only their ab- 
solute values. 

Figure QJa) shows the relative matrix elements 
Rrei (nSi/2 — > E'Pi/ 2 ) for bound-free transitions in ru- 
bidium atoms. The horizontal axis is the principal quan- 
tum number of the nS states, while the vertical axis 




FIG. 4. (Color online). Density plots of the relative matrix 



E'Pi/2 transitions; 



elements for (a) Rb bound- free nSi/2 

(b) arbitrary Rb ES1/2 — > E'Pi/ 2 transitions including dis- 
crete and continuum spectra; (c) Rb bound-bound nSi/2 — ► 
n'Pi/2 transitions; (d) Rb bound-free nP 1 / 2 —> E'Si/2 tran- 
sitions. 



is the binding energy E' of the continuum P states (in 
atomic units). From Fig.^a) one finds that for a bound- 
free transition between the Rydberg 40S , 1 / 2 state and 
the continuum E'Pi/2 state with E' = 0.18 the relative 



matrix element is R re i \nS\ 



/2 



E'P 1 



/2J 



0.02. Ac- 



cording to Eq.((T]) and Table 1 in the appendix A the 
energy of the 40Si/ 2 state is E n — —3.68 x 10~ 4 , and 
the energy difference is \E' — E n \ = 0.18. Then the ab- 
solute radial matrix element can be found from Eq. : 
R(40S 1/2 -> E'Pi/2) = 6.1 x 10" 4 . To highlight Cooper 
minima, the regions where the relative matrix element 
Rrei (nSi/2 — > E'P 1 / 2 ) falls below 0.02 are filled by black. 
The interesting feature of Fig. 01a) is the presence of the 
two sharp Cooper minima at E' = 0.09 and E' = 0.21 
for n > 10. 

Figure 0Jb) shows the relative matrix elements for 
all possible bound-bound, bound-free and free-free 
Rrei (S'1/2 Pi/ 2) transitions in rubidium atoms, plot- 
ted in the energy scale. The horizontal axis carries the 
binding energy of S states while the binding energy of 
P states is given in the vertical axis. For a bound-free 
transition between the QP1/2 state with E = —0.045 
and the continuum Si/2 state with E' = 0.15, the rel- 
ative matrix element is R re i (E n P\/2 E'Si/ 2 ) = 0.8. 
Then Eq. ([2]) gives the absolute radial matrix element 
R (E n P 1/2 -> E'Si/2) = 0.82 (the energy difference is 
\E' — E\ = 0.195). The same procedure can be applied 
to calculate radial matrix elements for all bound-bound, 
bound- free and free- free transitions. 

A prominent Cooper minimum is observed in 
Fig. BJa), (b) for the bound-free nSi/ 2 — > E'Pi/2 an d 
free-free ESi/ 2 — > F'P\/2 transitions in rubidium. For 
the nS\/2 states with E n > —0.03 (corresponding to 
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FIG. 5. (Color online). Density plots of the relative matrix 
elements in Rb atoms for (a) bound- free nPs/2 E'D^/ 2 
transitions; (b) all EP3/2 — > E'D 5 / 2 transitions; (c) bound- 
bound nP s /2 — ► n " D5/2 transitions; (d) bound-free nD 5 / 2 — ► 
E'P 3 /2 transitions; (e) bound-free nD 3 y 2 — ¥ E'F 5 / 2 transi- 
tions; (f) all ED3/2 — ¥ E'F 5 /2 transitions; (g) bound-bound 
nD 3 /2 —> n ' 'F5/2 transitions; (h) bound-free nF 5 / 2 E'D 3 / 2 
transitions. 



FIG. 6. (Color online). Density plots of the relative ma- 
trix elements in Cs atoms for (a) bound- free nSi/2 — ¥ E'P±/2 
transitions; (b) all ES1/2 — > E'Px/2 transitions; (c) bound- 
bound nS\/2 —¥ n'Px/2 transitions; (d) bound-free nP 1( / 2 — > 
E'Si/2 transitions; (e) bound-free nP^/2 — ¥ E'D 5 /2 transi- 
tions; (f) all -E-P3/2 — ¥ E'D 5 / 2 transitions; (g) bound-bound 
nP3/ 2 n'D 5 / 2 transitions; (h) bound-free nD 5 / 2 — > E'P 3 / 2 
transitions. 



n > 7) the relative matrix elements fall down below 0.02, 
while for the lower nSx/2 states the minimum is not so 
sharp. 

The relative matrix elements R re i {E n P\/ 2 En'Si/2) 
for bound-bound transitions in rubidium are presented in 
Fig-U^c). Since the relative matrix elements slowly vary 
with n and n', accurate calculation of the radial matrix 
elements from the data of Fig. @Ic) is possible as de- 
scribed earlier. For example, R re \ (27S , 1 / 2 — ¥ 8OP1/2) = 
0.45, E 27S = -8.776 x 10~ 4 , E 80P = -8.36 x 10~ 5 and 
the energy difference is \E n > - E n \ = 7.94 x 10~ 4 . From 
Eq. © one finds R (27S 1/2 -> 80P 1/2 ) = 0.342. 

Relative matrix elements R re i {nP\/ 2 — ^ E' Sx/2) for 
bound-free transitions in rubidium atoms are shown in 
Fig. @[d). As an example, the relative matrix element 



is R rel (60Pi /2 -> E'Sx/2) = 0.84 for E' = 0.2; then the 
radial matrix element is R (60P 1/2 -> E'S 1/2 ) = 0.012. 

Figure [5] displays the relative matrix elements 
Rrei (P5/2 -> P5/2) and R re i (P3/2 ~> P5/2) for rubid- 
ium atoms in the same way as in Fig. 01 Relative matrix 
elements for other fine-structure components of the ru- 
bidium P and D states are not presented, since the dif- 
ference between them is too small to be distinguishable 
on our density plots. Cooper minima are observed for 
bound-free nP 3 / 2 — ¥ E'D 5 / 2 transitions with E 1 s» 0.37 
[Fig. EJa)], free-free PP 3/2 -> E'D 5/2 and ED 3/2 P 5/2 
transitions [Fig. EJb)], and bound-free nF 5 / 2 —¥ E'D 3 / 2 
transitions with E' w 0.17 [Fig. EJh)] . 
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nD^ D^-state energy [a.u.] 



FIG. 7. (Color online). Density plots of the relative matrix 
elements in Cs atoms for (a) bound-free nD 3 / 2 — ► E'F^/2 
transitions; (b) all ED s / 2 E'F 5 / 2 transitions; (c) bound- 
bound nD :i / 2 —> n'F 5 /2 transitions; (d) bound-free nF 5 / 2 — ► 
E'D 3 / 2 transitions. 



B. Cesium 

Figure [6] shows relative matrix elements 
Rrei (5*1/2 -> Pa/2) and R re i (P 3/2 -> D 5/2 ) for cesium 
atoms. Relative matrix elements R re i {D3/2 ^ P5/2) 
for cesium atoms are presented in Fig. [7] Two Cooper 
minima are observed for bound- free nSi/ 2 — > E'P 1 / 2 
transitions with E' w 0.06 and 2?' « 0.17 in Fig. [BJa). 
The Cooper minima are also noticed for the bound- 
free nD 5 / 2 EP 3 / 2 transitions with E' ~ 0.28 
[Fig. IHh)], bound-free nF 5 / 2 — > E'D$m transitions with 
£?' = 0.19[Fig.[^h)], and bound-bound nD a/2 -> n'F 5/2 
transitions with n < 23 [Fig. Efg)]. 

Cooper minimum in the discrete spectrum was first 
discussed in Ref. [f|. In Ref. [IH a continuation of this 
minimum in the bound- free nD 3 / 2 — > E'F 5 / 2 transitions 
was found, with the energy 2?' ss 10 -3 being close to the 
ionization threshold. These features were reproduced in 
our calculations, but they are not shown due to the large 
energy scale of our density plots. The near-threshold 
Cooper minima for Rydberg states of alkali-metal atoms 
were also discussed in Ref. 111. 




nS S-state energy [a.u.] 





nP P-state energy [a.u.] 



FIG. 8. (Color online). Density plots of the relative ma- 
trix elements in Li atoms for (a) bound-free nS — > E'P 
transitions; (b) all ES — > E'P transitions; (c) bound-bound 
nS n'P transitions; (d) bound-free nP — > E'S transitions; 
(e) bound-free nP -» E'D transitions; (f) all EP ->■ E' D 
transitions; (g) bound-bound nP — > n'D transitions; (h) 
bound-free nD — > E'P transitions. 

transitions [Fig. [5fa)]. This minimum has been found 
earlier in Ref. [lOj in an experimental study of the far- 
infrared transitions between Rydberg states. Later on we 
have shown that such minimum can also appear in the 
BBR- induced transitions (50j . 



C. Lithium 

For lithium atoms the relative matrix elements ob- 
tained for R re i (S — > P) and R re i (P — > D) transitions 
are shown in Fig. [31 and for R re i (D — > F) transitions 
in Fig. 13 Fine structure is neglected due to the small 
fine splitting. One may see that the matrix elements of 
the bound-free S —> P transitions slowly decrease as the 
energy of the continuum state grows. A Cooper mini- 
mum is observed for the bound-bound nSi/ 2 — ¥ n'P\j 2 



D. Sodium 

For sodium atoms the relative matrix el- 
ements obtained for Rrei ipi/2 Pi/ 2) an( l 
Rrei (-P3/2 — * -O5/2) transitions are shown in Fig.[l0l and 
for Rrei (-D3/2 — > F5/2) transitions in Fig. [TTJ A Cooper 
minimum is observed for the bound-free nS\/ 2 — > E'Pi/2 
transitions at E' « 0.06 [Fig. ITOT a)]. bound-free 
nD 5/2 -> E'P 3/2 transitions at E' w 0.05 [Fig. HSfh)] , 
bound-free nP 3 / 2 —> E'D 5 / 2 transitions at E' w 0.25 
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FIG. 9. (Color online). Density plots of the relative ma- 
trix elements in Li atoms for (a) bound-free nD — > E'F 
transitions; (b) all ED — > E'F transitions; (c) bound-bound 
nD —¥ n' F transitions; (d) bound-free nF — > E'D transitions. 



[Fig. [He)] 
tions at E' 
and nD 5 /2 



and bound-free nFc 



5/2 



E'D 3 / 2 transi- 



0.29 [Fig. UTTd)]. For nS 1/2 E'P 1/2 
E'P 3 /2 transitions in sodium a Cooper 
minimum is found to be close to the ionization threshold. 



E. Potassium 



potassium atoms the relative matrix el- 
obtained for R re i (Si/2 —> Pi/ 2) an d 
{P3/2 -C5/2) transitions are shown in Fig. [TSJ 
for R re i (D 3 / 2 —> F5/2) transitions in Fig. [T3J 



For 
cments 

Rrei 

and 

Interesting features are observed in the radial matrix 
elements of bound-bound and bound- free transitions. 
The two minima have been located in the relative 
matrix elements of the bound- free nSi/2 — > E'P\/2 
transitions at E' « 0.05 and E' « 0.45 [Fig. HDJa)]. A 
Cooper minimum is also observed for nP 3 / 2 F'D 5 / 2 
transitions at E' « 0.29 [Fig. [T27g)]. A minimum in 
the discrete spectrum has been found [Fig. H2T c)] , which 
is similar to the lithium S — > P and cesium D — > F 
transitions. This minimum in potassium atoms was first 
discussed by Theodosiou [||. Finally, for the bound- free 
E'D 3 /2 transitions in potassium, a Cooper 



nF ; 



5/2 

minimum is registered at E 

pig. nid)]. 



0.18 and E' 



0.38 



IV. CONCLUSION 

The quasiclassical model developed by Dyachkov and 
Pankratov [33|, [H| can be used for fast and reliable cal- 
culations of the radial matrix elements for bound-bound, 
bound-free and free-free transitions between arbitrary 
states of alkali-metal atoms. We have demonstrated this 
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FIG. 10. (Color online). Density plots of the relative ma- 
trix elements in Na atoms for (a) bound-free nS\/2 — > E'Pi/ 2 
transitions; (b) all ES1/2 — ► E'Pi/2 transitions; (c) bound- 
bound nSx/2 —5- n'Pi/2 transitions; (d) bound-free nPi/2 — > 
E'Si/2 transitions; (e) bound-free nP^/2 —¥ E'D 5 / 2 transi- 
tions; (f) all EP 3 /2 E'D 5 /2 transitions; (g) bound-bound 
n P'i/2 n 'D 5 / 2 transitions; (h) bound-free nD 5 / 2 — ¥ E' P s 



3/2 



transitions. 



by perfoming the numerical calculations of the radial 
matrix elements for transitions between S, P, D and F 
states with the energies E < 0.5 (in atomic units) in 
all alkali-metal atoms. Our results on radial matrix el- 
ements are in good agreement with numerical calcula- 
tions in the Coulomb approximations [TTj | . Our theoret- 
ical results [lij are also consistent with the experimen- 
tal measurements of the effective lifetimes of Rydberg 
states [H, [2t| , oscillator strengths [22|, [23[ , and photoion- 
ization cross-sections Q. Our approach allowed us to 
reveal several unknown Cooper minima, both in the dis- 
crete and continuum spectra, which would be interesting 
to confirm experimentally. Reliability of the quasiclassi- 
cal model for study of the Cooper minima is verified by 
good agreement with calculations of Ref. @ for bound- 
bound transitions and satisfactory agreement with calcu- 
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FIG. 11. (Color online). Density plots of the relative matrix 
elements in Na atoms for (a) bound-free nD 3 / 2 ~ ► E'F B / 2 
transitions; (b) all ED 3 / 2 E'F s / 2 transitions; (c) bound- 
bound nD 3 / 2 — ► n'F B / 2 transitions; (d) bound-free nF 5 / 2 ~^ 
E D3/2 transitions; 



lations of Ref. [8J for bound-free transitions. We conclude 
that the quasiclassical model of Dyachkov and Pankra- 
tov is a universal method for systematic calculation of 
the radial matrix elements for transitions between ex- 
cited states of alkali-metal atoms. 
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FIG. 12. (Color online). Density plots of the relative ma- 
trix elements in K atoms for (a) bound-free nSi/ 2 — s- E' P 1 / 2 
transitions; (b) all ESi/ 2 — s> E 1 P^/ 2 transitions; (c) bound- 
bound nSi/2 —5- n'Pi/2 transitions; (d) bound-free nP 1 / 2 — > 
E'Si/2 transitions; (e) bound-free nP%/ 2 — » E'D$/ 2 transi- 
tions; (f) all EP- i j 2 — > E'D$/ 2 transitions; (g) bound-bound 
T1P3/2 — > n'D 5 / 2 transitions; (h) bound-free nD 5 / 2 — > E'P^/ 2 
transitions. 



Appendix A: The quantum defects. 



The values of the quantum defect can be obtained by 
fitting the experimentally measured energies to Eq.([T]) 



Ml (n) = a' L + + 



+ 



"eff 



L eff 



L eff 



"eff 



(Al) 



where a 



hi u Li L ' U L' 



are the Rydberg-Ritz fit- 
ting coefficients. The quantum defects can also be ex- 
pressed through the modified Rydberg-Ritz coefficients, 
which were tabulated for alkali-metal Rydberg atoms in 
Ref. (HI: 



[i L (n) = a L 



(n - a L ) 



(n - a L ) 



(n - a L ) (n - a L ) 



(A2) 



The difference between the fitting coefficients in 
Eqs. (|Al[) and (|A2I) is small and lies within the measure- 
ment uncertainty [5l| . The most recent experimental val- 
ues of the modified Rydberg-Ritz coefficients [57J avail- 
able for alkali-metal Rydberg atoms are listed in Table 1 
of Appendix. The data are taken from Refs. (5ll - [57| . 

We have compared our calculations with our recent 
experimental data for diffuse series in rubidium. Good 
agreement between experiment and theory is observed. 
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TABLE I. Quantum defects of alkali-metal Rydberg states 







Sl/2 


Pl/2 




D 3/2 


D 5/2 


F 5 /2 

F7/2 


Li m 


a 


0.39951183 


0.04716876 




0.00194211 




0.00030862 


fe 


0.02824560 


-0.02398188 




-0.00376875 




-0.00099057 


c 


0.02082123 


0.01548488 




-0.01563348 




-0.00739661 


d 


-0.09793152 


-0.16065777 




0.10335313 






e 


0.14782202 


0.33704280 










Ma 


a 


1.34796938(11) 


0.85544502(15) 


0.85462615(12) 


0.014909286(97) 


0.01492422(16) 


0.001453* 


b 


0.060989(16) 


0.112067(86) 


0.112344(67) 


-0.042506(35) 


-0.042585(43) 


0.017312* 


c 


0.019674(17) 


0.0479(13) 


0.0497(10) 


0.00840(31) 


0.00840(39) 


-0.7809* 


d 


-0.001045(354) 


0.0457(43) 


0.0406(34) 






7.021* 


K [H 


a 


2.1801985 


1.713892 


1.710848 


0.276970 


0.2771580 


0.010098 


b 


0.13558 


0.233294 


0.235437 


-1.024911 


-1.025635 


-0.100224 


c 


0.0759 


0.16137 


0.11551 


-0.709174 


-0.59201 


1.56334 


d 


0.117 


0.5345 


1.1015 


11.839 


10.0053 


-12.6851 


e 


-0.206 


-0.234 


-2.0356 


-26.689 


-19.0244 




Rb [HHH 


a 


3.1311804(10) 


2.6548849(10) 


2.6416737(10) 


1.34809171(40) 


1.34646572(30) 


0.0165192(9) 
(F 6/2 ) 














0.0165437(7) 

(F7/2) 




0.1784(6) 


0.2900(6) 


0.2950(7) 


-0.60286(26) 


-0.59600(18) 


-0.085(9) (F 5/2 ) 


b 












-0.086(7) (F 7/2 ) 


cs HH 


a 


4.04935665(38) 


3.59158950(58) 


3.5589599** 


2.4754562** 


2.46631524(63) 


0.03341424(96) 


b 


0.2377037 


0.360926 


0.392469** 


0.009320** 


0.013577 


-0.198674 


c 


0.255401 


0.41905 


-0.67431** 


-0.43498** 


-0.37457 


0.28953 


d 


0.00378 


0.64388 


22.3531** 


-0.76358** 


-2.1867 


-0.2601 


e 


0.25486 


1.45035 


-92.289** 


-18.0061** 


-1.5532 





* 51' 

**lj5' 




D -state energy [a.u.] 



FIG. 13. (Color online). Density plots of the relative matrix 
elements in K atoms for (a) bound-free nD 3 / 2 — ► E'F 5 / 2 tran- 
sitions; (b) all ED3/2 — > E'F 5 /2 transitions; (c) bound-bound 



nD 



3/2 



% ' 'F5/2 transitions; (d) bound-free nF 5 / 2 ~ > E'D £ 



3/2 



transitions. 



For the calculation of the dipole matrix elements 
of bound- free and free-free transitions, the expressions 
for the quantum defects of Rydberg states, given by 



Eqs.(|AT 
tinuum 



and (|A2[) . must be extrapolated to the con- 



li L (E) = a L + b L x (-2E) + 
+c L x (-2E) 2 +d L x (-2E) S 



(A3) 



Here E is energy of the continuum state. We have found 
that calculations of the radial matrix elements for bound- 
free and free-free transitions are sensitive to the way of 
extrapolation of the quantum defects into continuum, 
especially in the regions of the Cooper minima. This 
was also discussed earlier in Ref. [49|]. The most re- 
cent data for the rubidium quantum defects [54l [55j con- 
tain only two coefficients cll and b^, while for the other 
alkali- metal atoms up to five terms of Eq. (|A2[) have been 
published [5ll - l53l . [5a . [57| . However, our test for using 
the higher-order polynomial approximations gave incor- 
rect values of the quantum defects at large energies of 
the continuum states, due to the well known Runge's 
phenomenon [59j of oscillation of the interpolation func- 
tion near the edge of the interpolation region, when the 
higher-order polynomial approximation is used. There- 
fore, in our calculations we applied a linear extrapolation 
of quantum defects to the continuum with only the first 
two terms of Eq. (|A"3|) . 
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Appendix B: Dyachkov-Pankratov quasiclassical 
model 



A quasiclassical model of Dyachkov and Pankratov was published in their original papers [3J, [3J] . Here we sum- 
marize the main formulas, which are used to calculate the relative matrix elements for bound-bound, bound-free and 
free-free transitions between \E,L) and \E',L') states and transition frequency to = E' — E, wherei?' > E. The 
difference of the quantum defects Afi = fiL' (E') — (E). 

The Keppler motion of quasiclassical electron is determined by the mean energy: 



E r 



E + E' 



(Bl) 



The cases of E c < (finite mean orbit) and E c > (infinite mean orbit) must be considered separately. If E c < 
(bound-bound transitions and bound-free transitions to the continuum states with E' < \E\), the parameters 7*, 7 
and mean quantum number v c are defined as: 



\E + E'\ 3/2 '' 



7 = int [7* + A/J. + 0.5] - A/z, 
vc = (7/ w ) 1/3 ■ 



(B2) 



Here wit[a;] means integer part of x. In the quasiclassical model the dynamics of the electron is defined by the 
arithmetic mean orbital momentum l c and eccentricity of the mean elliptic orbit of the electron e c : 



L 



L + L' + l 



2 \ 1/2 



(B3) 



The relative matrix element is then expressed as: 

- - vW 3 



Rrel {EL — > E L 



0.4108 



[U~ (£7) cos (vrA/i) — Vj (£7) sin (irAfi)] 



f/ 7 = J' (£7) + A/ — J 7 (e 7 ) 



V C E 



V 1 = EL (£7) + Ale- 



ves 



E 1 ( £7 ) 

7T7 



1 -£ 



(B4) 



Here J 7 and E 1 are the Anger and Weber functions, respectively, and E' are their derivatives with respect to 
argument. 

If E c > (bound- free transitions with E' > \E\ and free-free transitions), the relative matrix element is given by: 

2 2/3 

R rel (EL -> E'U) = [Pi (£7) cos (ttA^) - Q 7 (£7) sin (vrA^)] , 



P 1 = -g' (£7) + AZ-^-«? 7 (£7) , 



Q 7 = h' (£7) + Al- 



lies 



hy (£7) 



7T7 



(B5) 



Here r) c = l/y/2E~ c , 7 = f^w and £ = (l + l 2 /Ve) 1/2 - 

Functions g 7 and /i 7 are expressed through the Hankel Anger J and modified Bessel / functions: 



foy (y) = 



1 



sinh (ttj) 



Ji~t (-%) - ^ ex P 07/2) [i»J/ (7) + (7)] 



(B6) 
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In the case of E c ~ an asymptotic expression for relative matrix element can be used: 

41/3 



Rrel {EL E L 



{[S (x) + /3S 1 (x)] cos (ttA^) - [T (x) + /3Tj (x)] sin (ttA^)} 



(B7) 



0.4108 

Here x = 3 , /3 = 2E c (2/uj) 2 ^ 3 , and functions 5(x), Si (x), T(x), 7\ (x) are expressed through the Airy 

Ai (x) and Bi (x) functions, their derivatives and hypergeometric function \F 2 : 



S (x) = Al ■ x 1/2 ■ Ai (x) - Ai' (x) , 
T (x) = Al • x 1/2 • Gi (x) - Gi' (x) + 

1 (x) = (l - 6AI • x 3/2 + 4x 3 ) Ai ( 



x 

2tt' 



S 



Ti (x) = — - 6AZ • x 3/2 + 4x 3 J Gi (x) 

Gi (x) = -Bi (x) - — X F 2 f 1; -, -; — 
V 7 3 W 2tt I ' 3' 3' 3 



(l - AZ • x 3 / 2 ) Ai' (x) , 

f 1 - Al- x 3 / 2 ) Gi' (x) + J-Al ■ x 1 ' 2 

V / 20?T 



21 



-Bi (x) x 
— ] 

7T 



4 5 x 3 
3' 3' "9 



3x^ 
40tt j 



Gi' (x) = — ^ - Z1 F 2 ( 1; -,-;-) - — 1F 2 ( 2; -,-;-] — 



x 

2%' 



(B8) 



We note that in the original paper [331 ] there were misprints in the last two terms in the expression for T\ (x). Here 
we present the corrected formula kindly provided by L. G. Dyachkov [60( ] . 

To simplify numerical calculations with generic mathematical codes and software, functions J 7 , F 7 , g 7 ,/i 7 and 
their derivatives can be expressed via more commonly used reg ularized hypergeometric functions 1-F 2 , hypergeometric 
functions 1-F 2 , Hankel H and modified Bessel I functions [611 : 



i; ^(3-7)^(3 + 7); ~ 



/7T7 



COS ^— ) !^ 



7 7 z 

l; i + tt; — 



21 



Ry{z) 

+ 



ii 5(3-7), 5(3 + 7); 



7 7 z 

1;1 -i' 1 + i ; -T 



2; 2- -. 2 ' 

2' 



3—7 3+7 



2' ~4 

2 



Z sinf^ ,F 2 



A 2 / 



1: 1 



2: 2-?, 2 + ?: -- 



1 



5 ~ 

1 



3-7 3+7 
' 2 ' 2 : 



1: 1 



3-7 



3-7 



3 + 7 _z 2 
2 5 4 



3 + 7 _ z 2 
2 ' 4 



h y ( z ) = cosh(7T7) |-^exp (Jy) [I-i-h ( z ) + h-vy 0) + /-i+ i7 (2) +7i +i7 (z)]| + 

2 v 2 / v 2 2 4 j 



-z coshf^F 2 f 2 ;2-^, 2 + !2 ; f: 
2 V 2 / V 2 ' 2 4 



n F 2 1:1 



3 — 27 3 + 17 z' 



(B9) 
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